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Abstract 

We study the dynamical symmetries of a class of two-dimensional superintegrable sys- 
tems on a 2-sphere, obtained by a procedure based on the Marsden-Weinstein reduction, 
by considering its shape-invariant intertwining operators. These are obtained by gener- 
alizing the techniques of factorization of one-dimensional systems. We firstly obtain a 
pair of noncommuting Lie algebras su(2) that originate the algebra so(4). By considering 
three spherical coordinate systems we get the algebra u(3) that can be enlarged by 're- 
flexions' to so(6). The bounded eigenstates of the Hamiltonian hierarchies are associated 
to the irreducible unitary representations of these dynamical algebras. 



1 Introduction 

It is well known that integrable Hamiltonian systems play a fundamental role in the descrip- 
tion of physical systems, because their many interesting properties both from mathematical 
and physical points of view. Many of these Hamiltonian systems have proved to be of an ex- 
traordinary physical interest e.g. harmonic oscilator, Kepler problem, Morse pQ, Posch-Teller 
[2], Smorodinski-Winternitz 011] > Calogero [5] and Sutherland jH] potentials. 

In this paper we present a class of integrable Hamiltonian systems that allow us to gener- 
alize the intertwining transformations for one-dimensional (ID) systems (7j to higher dimen- 
sions. 

These Hamiltonian systems are superintegrable, i.e., they have more than N constants 
of motion, being N the dimension of the configuration space for the Hamiltonian system. 
Although these motion integrals are not all of them in involution they determine more than 
one subset of N constants (in all the cases one of them is the Hamiltonian) in involution. The 
system is said to be superintegrable in the sense of [Hj or maximally superintegrable if there 
exist 2 N — 1 invariants well defined in phase-space. 

Using the Marsden-Weinstein reduction procedure [S] we construct such classical systems 
starting from a free Hamiltonian lying in an N-dimensional homogeneous space of a suitable 
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Lie group, whose action allows to us to calculate a momentum map that assure the inte- 
grability, or even superintegrability of the reduced systems (!]. From an opposite point of 
view, there are good reasons to suspect that any integrable system may be constructed as a 
reduction from a free one |1U| . 

For the corresponding quantum systems we present here a generalization to higher di- 
mensional spaces of the intertwining transformations for ID factorizable systems. These ID 
systems have dynamical Lie algebras of rank one generated by the intertwining operators 
0. By using a concrete superintegrable Hamiltonian system with underlying symmetry the 
Lie algebra su(3) we find that its dynamical symmetry can bee enlarged to so(6). However, 
these results can be implemented to higher dimensional systems of the same class, and also 
can be helpful in the study of other kinds of integrable systems using algebraic methods 

[mcaiiaimEiiiiii. 

In sections |2j 01 and 0] we will introduce a two-dimensional superintegrable system and 
find some separable solutions by standard procedures (Hamilton- Jacobi equation). Next, in 
sections and we study the corresponding Schrodinger equation which is factorizable in two 
ID equations. We construct some sets of intertwining operators closing the Lie algebras u(3) 
and so(6) by taking also into account different separable coordinate systems as well as discrete 
symmetries. We characterize the eigenfunctions of the Hamiltonian hierachies, obtained from 
the intertwining operators, as irreducible unitary representations of these dynamical algebras. 

2 Superintegrable SU (p, g)-Hamiltonian systems 

Let us consider a free Hamiltonian H = 4 g^ v p^p u (f-i, v = , • • • , n = p + q — 1 and the bar 

stands for complex conjugate) defined in the configuration space su(p~\ P q)y.U{i) • This is an 
hermitian hyperbolic space with metric and coordinates y^ 1 " G C, verifying g^ u y^ y v = 1 
(by Pfj, we denote the conjugate momenta). The geometry and properties of this kind of 
spaces are described in ^1 ^3] . 

Using a maximal abelian subalgebra (MASA) of the Lie algebra su(p, q) ^S] the reduction 
procedure allows us to obtain a reduced Hamiltonian, H = \g^ v p s v p s » + V(s), lying in 
the corresponding reduced space (a homogeneous SO(p, g)-space), where V(s) is a potential 
depending on the real coordinates satisfying the constraint g^ s u = 1. 

The set of complex coordinates y^ is transformed by the reduction into a set of ignorable 
variables (which are the parameters of the transformation associated to the MASA of 
u(p,q) used in the reduction) and the actual real coordinates 

If , fj, = 0, • • • , n, is a basis of the considered MASA of u(p, q), formed by pure imaginary 
matrices (this is a basic hypotesis in our reduction procedure), the relation between old (y M ) 
and new coordinates (x M , s^) is 

y» = B(x)t s v , B(x) = exp (xP . 

This relationship assures the ignorability of the x coordinates (the vector fields corresponding 
to the MASA are straightened out in these coordinates). The Jacobian matrix, J, of the 
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coordinate transformation is given explicitly by 



j_ a^y) _ (A B\ djf_ 

J ~ d(x,s) ~ \A Bj ' Av ~ dx» ~ [ v) f> V ■ 

The Hamiltonian calculated in the new coordinates is written as 

H = c(~g^ PflPv + V(s)\ , V{s) = p T x {A^ KA)- l Px , 

where p x are the constant momenta associated to the ignorable coordinates x and K is the 
matrix defined by the metric g. A detailed exposition of this construction procedure of this 
family of superintegrable systems can be found in [20] . 



3 A classical superintegrable w(3)-Hamiltonian system 

To obtain the classical superintegrable Hamiltonian associated to the unitary Lie algebra 
su(3) using the reduction procedure sketched in the previous section we are going to proceed 
in the following way |21| . 

Let us consider the basis of su(3) determined by 3 x 3 matrices X\ , ■ ■ ■ ,X&, whose explicit 
form, when using the metric K = diag(l, 1,1), is 

li 0\ /0 \ / 1 0\ /0 i 0\ 

X x = -i X 2 = i X 3 = 1-1 X A = li 
\0 0/ \0 -i) \ 0/ \0 0/ 



/ 1\ /0 A /0 0\ /0 0\ 

X 5 = X 6 = X 7 = 1 X 8 = i . 
\-l 0/ \i 0/ \0 -1 0/ \0 i 0/ 

In the compact case, as here with su(3), there is only one MAS A. This is the Cartan 
subalgebra, generated by the matrices diag(z,— i,Q) and diag(0, i, — i). However, we shall 
work, in order to facilitate the computations, with the algebra u(3) instead of su(3). Hence, 
we shall use the following basis for the corresponding MASA in u(3) 

y = diag(i,0,0), Fi = diag(0,i,0), Y 2 = diag(0, 0, i). 

The actual real coordinates s are related to the complex coordinates y by = s M e 1 x ^ (/i = 
0, 1, 2). The Hamiltonian can be written as 

H=U P 2 o+ P i+p 2 2 )+V(s), V(s)=4 + 4 + 4> (3-1) 
Z Sq s± s 2 

lying in the 2-sphere (s ) 2 + (si) 2 + (s 2 ) 2 = 1. 

To see that the system, so obtained, is superintegrable it is necessary to construct its 
invariants of motion. In this case we obtain three invariants 

Rav = {suPv - s u pu) 2 + ( m»— + m u — ) , fj,< v, n = 0, 1, u = 1,2. (3.2) 

V Sf, s v ) 
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Note that only two of them are in involution at the same time (being one of them the 
Hamiltonian) , so the system is sup erintegr able in the sense of [S]. The sum of these invariants 
is the Hamiltonian up to an additive constant. 



4 The Hamilton-Jacobi equation for the w(3)-system 

The solutions of the motion problem for this system can be obtained solving the corresponding 
Hamilton-Jacobi (HJ) equation in an appropriate coordinate system, such that the Hamilton- 
Jacobi equation separates into a system of ordinary differential equations. 

The 2-sphere can be parametrized on spherical coordinates (<j>i,(f>2) around the S2 axis 

so = cos02 cos 0i, si = cos02 sin0i, S2 = sin 02 (4-1) 
where 0i G [0, 2ir) and 02 £ [tt/2, 3ir/2]. Then, the Hamiltonian (|3.1[) is written as 



cos 2 02 \ cos 2 0i sin 2 0i / sin 2 02 

The potential is periodic and has singularities along the coordinate lines 0i = 0, vr/2, tt, 37r/2, 
and 02 = vr/2, 3tt/2, and thre is a unique minimum inside each regularity domain. 

The invariants ()3.2j) (denoted by /) can be written in terms of the basis {Xi , • • • , Xg} as 
h = Xl + X\ , I 2 = Xl + Xf and I 3 = X% + X%. The u(3)-Casimir is 

C = iXf + 2{Xi , X 2 } + 4X 2 + 3/i + 3/ 2 + 3/ 3 , 

where the three first terms are the second order operators in the enveloping algebra of the 
compact Cartan subalgebra of u(3). 

The Hamiltonian is rewritten as H = I\ + I2 + ^3 + constant, where by Ii we denote the 
invariant I{ but rewritten in spherical coordinates. So, we have 

1 2 m 2 mi2 



2 ripi cos 2 0i sin 2 0i ' 

- / ^ 2 ■ 2 1 m o \ 1 2 1 ( 1 2 m % 



tan 2 y-p^ sin 0i + — ^-j + cos 0i ^-p^ + 

+ i^P<t>iP<j>2 sin 2 0i tan 02, 
2 1 ( 1 2 2 i m i \ • 2 i /" 1 2 m ^ 



'2 



tan 02 -p^ cos 0i + , 2 + sin 0i \-p (j)2 + 
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2'" 1 — Ti ' sin 2 0i ; ' " ri ' tan 2 02 

p^Pfc sin 2 0i tan 2 . 
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Now, the HJ equation takes the form 




+ 



+ 



cos z 



sin 2 



E. 



It separates into two ordinary differential equations taking into account that the solution of 
the HJ equation can be written as S((f>i, <fo) = 5i(^i) + £2(^2) — Et. Thus, 



dSi 



+ 



+ 



COS 2 i 



+ 



sin 2 (f>i 



sin 2 4>2 



+ 



a i 



COS^ 



OCX, 



OL2, 



where 02 = E and cti are the separation constants (which are positive). Each one of these 
two equations have the same form of those corresponding to the ID problem |21| . 

The solutions of both HJ equations are easily computed and can be found as particular 
cases in Ref. j!8J . A detailed analysis of them shows that all the orbits in a neighborhood 
of a critical point (center) are closed and thus, the corresponding trajectories are periodic 
(a direct consequence of the correspondence between extrema of the potential and critical 
points of the phase space). 

The explicit solutions, when we restrict us to the domain < <f>i,(f>2 < 7T /2, are 



COS 2 (j)2 



cos 



1 

2E 
1 

2q^ 



b 2 + 
61 + 



b% - AaiEcos2V2Et 



cos^ 



b\ - 4a\E 



1/2 



+ 2 v / aT[(&2 - E COS 2 4>2) COS 2 (f>2 



((62 cos (f>2 — 2ai) sin 2^j2a>i(3i 
1/2 cos 2V2ai7?i 



where b\ = ol\ + m^ — m\ and &2 = E + a\ — m^. 

Note that in this domain for the variables (4>i,4>2) the minimum for the potential corre- 
sponds to the point (^i = arctan mi/mo, 02 = arctan ^2/^0 + mi)). The value of the 
potential at this point is (too + m-i + 7712) 2 . Hence, the energy E is bounded from below, i.e. 
E > (mo + mi + m.2) 2 . 

As we mentioned before these results reflect essentially the (ID) su(2)-case. In fact, all 
systems we can construct using Cartan subalgebras can be described in a unified way as it 
was shown in [TH1 I2771I2T] . 



5 A quantum superintegrable w(3)-Hamiltonian system 



From the quantum point of view the Hamiltonian (j.S.lj) takes the form 

/r2 T 2 t 2\ Jo" 1 / 4 *1-V 4 

H = -(4 + J\ 2 + J 2 ) + ^—i- + + - 



{sof 



(5.1) 
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where (lo,hih) S M 3 , and Jj = —CijkSjdk- 

Also in the quantum case, the eigenvalue problem H & = E after substituting the 
coordinates l|4.1j) . takes the form of a separable differential equation 



-<9? + tan(0 2 )d0 2 + 



sin 2 (0 2 ) ' cos 2 (0 2 ) 



+ 



2 ig-1/4 g-1/4 
^ ^ cos 2 (0i ) ^ sin 2 (01 ) 



$ = (5.2) 



Taking the solutions separated in the variables <f>\ and 2 as <&(0i,0 2 ) = f((pi)g(4>2) after 
replacing in (|5.2j) we get the equations 



& I g ~ 1/4 g ~ 1/4 
* cos 2 (0i) sin 2 (0i) 

-9l +tan(0 2 )5^ 2 + 



f{<h) = <xf(<h), 
q - 1/4 



+ 



cos 2 (0 2 ) ' sin 2 (0 2 ) 



(5.3) 
(5.4) 



where a is a separating constant. 

These two (one- variable dependent) equations can be solved using the standard factoriza- 
tions obtaining polynomial solutions. Notice that the results obtained for the first equation 
will match in a certain way with those of the second one giving rise to degenerate levels. 



5.1 The factorization of the 0i-equation 

The ID Hamiltonian corresponding to the equation ()5.3j) in the variable <p\ can be factorized 
using the theory of factorizations by Infeld and Hull [7j . 

The second order differential operator at the l.h.s. of equation (|5.3f) can be written as a 
product of first order operators 

H (o) = A o A o + X o, 

where = ±8^ - (l + 1/2) tan^i + (h + 1/2) cot 0i, and A = (k +h + l) 2 . Also it is 
possible to construct a family of operators {^4i, A~ , A m , H? 1 ,}, m G Z, where 

A± = ±8^ -(l + m + 1/2) tan <f> x + (h + m + 1/2) cot <j) X , (5.5) 
A m = (I + !i + 2m + l) 2 , 



01 2 (l + m) 2 - 1/4 (f 1+ m) 2 -l/4 



(5.6) 



^ COS^^l SUV01 

Hence, we obtain a ID Hamiltonian hierarchy Q5.6|) . whose first element is Hty., satisfies 

^(m) = ^™ ^™ + A m = ^4 m _i ^m-l + -Vn-1- (5-7) 

From it we see that the operators A^ are shape invariant intertwining operators, i.e. 



^m#"(m) — H^ n+l ,A m , — H,KA^. (5i 
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Formally, the operators A^ acting on a Hamiltonian eigenfunction give another eigenfunction 
of a consecutive Hamiltonian in the hierarchy with the same eigenvalue, i.e. 

a- . __. -oy'/'i 4+ . n-/^>i __. <xy</>i 

where Tim is the eigenfunction space of H$\ . 

The discrete spectrum and the physical eigenstates of may be obtained, in principle, 
from the fundamental states f9* (and their eigenvalues) of all the Hamiltonians of the hier- 
archy {Ht\}. The fundamental states are determined by the equation A^f^ = whose 
solutions, up to a normalization constant, are 

= co B Zo+m+1 /2(0 1 ) s in^ +m+1 / 2 (^) , (5.9) 

with eigenvalues A m = (Iq + /i + 2m + l) 2 . 

The excited eigenfunction fJjK of H$\ can be obtained from the ground eigenstate f? * 
of Hjpy (both with the same eigenvalue) applying consecutive operators A + 

f m — 4 + 4 + ... 4+ f° 

(5.10) 

= N sin /l+1 / 2 (0i) cos' 0+1 / 2 (0i ) Pr A,) [cos(2 </»!)], 

where Pn a ' b \x) are Jacobi polynomials and N a normalization constant. The spectrum of 
the Hamiltonian H£k (j5.3|) is given by 

a = A m = (I + I 1 + 2m + lf, m £ Z-° . (5.11) 



5.2 The dynamical algebras associated to the 01-factorization 

The shape invariant intertwining operators for the ID Hamiltonian hierarchy {Ht\} deter- 
mine some Lie algebras that we can characterize as follows. 

Let us define free-index operators A* 1 acting inside the total space ®mT~im from the above 

Af m) by 123123 

^4 + /(m+l) := 2^m/(m+l) f(m) 

A~f( m ) '■= ^ A mf(m) «C /( m+ i) (5.12) 

^/(m) := -|(fo+Zl+2m)/( m ) oc / (m) 
where /( m ) (or fr m )) denotes an eigenfunction of Ht m y We can rewrite Q5.7JI as 

[A, A*] = ±A ± , [A-,A + ] = -2A (5.13) 

assuming that the action is on any fr m y These commutators determine a Lie algebra su(2) 
whose Casimir element is C = A + A~ + A{A — 1). 

It can be proved (for more details see Ref. [2H1) that the eigenstates of the Hamiltonians 
H?m) Q5.9JI can be characterized, if Iq and l\ are positive o zero integer numbers, in terms of 
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the vectors of the irreducible unitary representations (IUR) of su(2), labeled by the parameter 
j such that 2j 6 Z— °. Effectively, the ground states f? ^ of Ht\ are characterized by 

A-f {m) = 0, Af (m) = -[(Z + /i + 2m)/2]/ ( ° m) . (5.14) 
Then, we identify (up to a normalization constant) 

f(m) = \jm, -jm), jm = {k + h + 2m)/2, 

where |j s) denotes the vectors of the IUR l j' of su(2). 

The excited states of Ht\ are obtained using expression (|5. lUf) . So, the eigenstate of the 
fc-th excited level of Htk is 

/(O) = lift + fc, -Jfc + fc) , 3k = (lo + h + 2k)/2, k = 0,1,2... 

Moreover, Htk (as well as any Ht\) can be expressed in terms of the su(2)-Casimir C acting 

on such representations by H?l = 4(C + 1/4). Hence, the eigenvalue equation for any of the 
excited states can be written as follows 

Hfcffc = 4(C + 1/4) |jo + k, -j + k} = 4(io + k + l/2) 2 |j + fc, -j + fc) 

= (Z + /i + 2fc + l) 2 / ( o), fc = 0,1,2,... 

However, there is an ambiguity due that different fundamental states ()5.9(l with values of 
Zo and l\ giving the same jo = (Zo + Zi)/2 would lead to the same jo-representation of su(2). 

Adding the diagonal operator D, Dfr m \ := (Zo — h)f( m ), to the generators of su(2) (|5.12|) 
we obtain u(2). The eigenstates of the Hamiltonian hierarchy are now completly characterized 
by the IUR's of u{2). However, different ii(2)-IUR's may give rise to (different) states with 
the same energy. 

On the other hand, the states of the Hamiltonian hierarchies, when Zo or 1% are not in 
Z-°, correspond to non-unitary representations of u(2), although they and their spectra are 
also given by formulae (|5.9|) and l)5.11j) . respectively. 

From a classification point of view it will be interesting to construct a new Lie algebra, 
obviously containing the subalgebra chain su(3) C u(3), such that only one of its IUR's 
characterizes all the eigenstates in the hierarchy with same energy. 

In order to build a such dynamical algebra we need to introduce a a two-subindex notation, 
in terms of the parameters (Zo, Zi), for the intertwining operators. The Hamiltonians H^\^ 

(j5.6j) will be denoted by H'fr its eigenfunctions by f(i ^), and the factor operators ^ in 
(j5.5j) will be rewritten as ; n. In this way, relation (J5.8|) can be expressed as 

A- TT<t>l __ ZT<Pl A- A+ JJ<t>l __ tt4>1 A + 

(loA) (lo,h) -"(Zo+l.ii+l) ^(Wi) ' (h,h) (lo+hh+l) (k,h) (Wi)' 



We can also define the free-subindex operators A^ 1 , A, D as in (|5.12j) . 

The fact that each two-parameter Hamiltonian Ht* ^ is invariant under the reflections 
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originates a second factorization ([27J |2H1 EH) via conjugation of the operators A^,A, D, 
I Q A±I Q = I± IqAIq = A, I DI = D, 

I 1 A ± I 1 = i=F, hAI x = -I, hDh = —D. 

The explicit form of the new operators is 

Af hM = ±0 01 + (io-1/2) tan^i + (h + 1/2) cot&, = --(-lo+h). (5.15) 

These operators {y^A 1 * 1 } close a second Lie algebra su(2), denoted su(2), which commutes 
with the previous su(2). Since, moreover, D and D essentially coincide with A and A, 
respectively, the new dynamical algebra is su(2) ®su(2) ~ so(4). 

The action of the so(4)-generators on a Hamiltonian ^ originates a 2D parameter 
so(4)-hierarchy {Hi _ n+m i 1+n+m } , to, n G Z, fixed by the initial values (Zo, ii) . Each energy 
level of this Hamiltonian hierarchy is degenerated and the eigenstates are characterized by 
so(4)-representations. 



5.3 The factorization of the 2 -equation 

The second equation (|5.4|) can be also factorized provided that the separation constant a is 
substituted by the eigenvalues obtained from the ^i-factorization a = X m = (Iq + l\ + 2m) 2 . 
The Hamiltonian associated to this (^-equation (|5.4|) is 

rth »2 uw | (/q + /i + 2to + 1) 2 II - 1/4 

^(0) = ~ d 4>2 + tan 2 )^ 2 + + . 2 . (5.16) 

w v cos z ((p2j sur((/) 2 ) 

It can be factorized in terms of two first-order differential operators as Ht^, = M^Mq + hq. 

This Hamiltonian fffe is the first element of the Hamiltonian hierarchy {ff^}, n E Z-°, in 
the 02-variable, whose elements can be written as 

H (n) = M nM~ +fJLn = M^_ X M+_ X + /i n _ 1; 

M± = ±d H - (fa + h + 2(to + 1) + n) tan(0 2 ) + (Z 2 + n + 1/2) cot(</> 2 ), 

/in = (h + h + h + 2n + 2m + 3/2)(Z 2 + /i + Z + 2n + 2m + 5/2). (5.17) 

The energy values E n are given by the factorization constant fj, n , i.e. E = fi n . The ground 
states g9 \ for this hierarchy are 

IX \ t\t I x \ii+in+2m+l • (± \l 2 +n+l/2 

9{n)VPV = N cos K<p2) sm(0 2 ) 2T T ' . 

The eigenfunctions gj^ of the initial Hamiltonian ijfe (|5.16() can be written as 

^ o) (0 2 ) = cos(0 2 ) /l+ ' o+2m+1 sin(0 2 )' 2+1 / 2 P^ 2+1 / 2 ' /l+ ' o+2m+1 )[cos20 2 ]. (5.18) 
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The index- free operators M^, denned in a similar way as A ± in (|5.12jl . close again a Lie 
algebra su(2). The eigenfunctions Q5.18JI are square-integrable, but the su(2)-representations 
are IUR provided that the parameters Iq, li,fo belong to Z-°. 

The new factorization leads to a degeneration of the energy levels indicating that the 
underlying dynamical symmetry could be larger than so(4). 

Finally joining both factorizations we obtain the square-integrable eigenfunctions sepa- 
rated in the variables (4>i, 02 ) of the Hamiltonian (|5.2j) 

where f^(<fii and g^(<f>2) are given by the expressions ()5.1U|> and l|5.18|) . respectively. Their 
corresponding eigenvalues E n ^ m (|5.17|) are degenerated for the values of m and n whose sum 
m + n keeps constant j2£] (see Figures ^ and |2| . 

6 Dynamical symmetries of the w(3)-Hamiltonian hierarchy 

The spectrum of the u(3)-Hamiltonian system (|5.1[) suggestes a bigger dynamical algebra 
of the Hamiltonian hierarchy. By introducing three sets of intertwining operators closing 
an algebra u(3) and using reflexion operators this algebra is enlarged to so(6). These three 
sets of operators are related with three set of spherical coordinates that we can take in the 
2-sphere submerged in a 3D ambient space with cartesian axes {so, s\, s?}. Since the axes 
(so, si, S2) play a symmetric role in the Hamiltonian (|5.1j) . we will take their cyclic rotations 
to get two other intertwining sets. These two set of spherical coordinates also separate the 
Hamiltonian 1)5.1(1 . 




Figure 1: Plot of the two IUR's so(4)-hierarchies, where each point represents a Hamiltonian. 
At the left it is the integer (Iq = 0, l\ = 0) and at the right the half-odd (Iq = 0, Zi = 1) 
hierarchy. The arrows stand for the intertwining operators. 
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Figure 2: The points of the plots represent eigenstates of the underlying half-odd Hamiltonian 
so(4)-hierarchy. At the left, the 4 eigenstates of the 5 ® 2 so(4)-IUR that share the energy 
E = (1 + l) 2 , and at the right the 16 eigenstates of the § (g> § so(4)-IUR with £ = (3 + l) 2 . 



6.1 The ti(3)-Hamiltonian hierarchies 

The Hamiltonian as well as the intertwining operators will be labelled now by three param- 
eters l±, I2) instead of two ones that we made used previously in section I5?2l 

The first set of sphere coordinates that we will use is that of coordinates (<f>i,<j>2) such 
that the coordinates {so> s i! s 2} are given by (|4.1|) . The Hamiltonian (|5.1[) characterized by 
the parameters t = (lo,h,h) will be now denoted by Hn j lt i 2 \, and the operators (|5.5|) will 
be rewritten with a three-fold subindex 

Af loMM) = ±8^ - (l + 1/2) tan 0i + (h + 1/2) cot^. 

The differential operators Q5.5|) depend only on the variable (pi, hence they do not act on the 
second separating variable 4>2- So, we have the intertwining relations 

A (l ,h,h) H Qo,h,h) = H (.lo+hh+l,h) A (l ,h,hV A (l ,h,h) H Qo+l,h+l,h) = H {hM2) A th,h,hY 
The global operators acting on eigenfunctions of these Hn Q j u i 2 \ defined as in (|5.12jl . 

A+ ®(lo+l,h+l,l2) := 2 A {l ,h,h)®(lo+i,h+i,h) K *(io,ii,i 2 )5 
A ~®Qo,h,h) '■= 2 A Ih->,ii,h)®( l o,hM K *(Jo+i,/i+i,i 2 )> 

A ®(lo,h,h) '■= --2~( lo+ll ^( l o,h,h)i 

close an algebra su(2) with commutators like ()5.13j) . Notice that now these operators are 
acting on the total wavefunction of the complete Hamiltonians like Hn^^^, not just on a 
factor function in one-variable. 

We will take the spherical coordinates choosing as 'third axis' s\ instead of S2, 
s 2 = cos(£ 2 ) cos(Ci), so = cos(£ 2 ) sin(£x), s x = sin(£ 2 ) • 



11 



The corresponding intertwining operators ^ ^ are denned in a similar way to ^ £ >. 
The explicit expressions for the new set in terms of the initial coordinates (|4.1[) are 

B Uoh h) = =l= ( sm( / , i tan 02001 + cos(/>iC^ 2 ) - (Z 2 +l/2) cos(/>icot0 2 + (l +l/2) sec^itan</> 2 . 
Their intertwining action on the Hamiltonians is 

B (lo,h,h) H ( l oM2) = H (lo+l,h,h+l) B (lo,li,h)> B (l ,h,h) H Oo+lM2+l) = H {k,h,h) B {l Q ,h,hY 

The 'global' operators, defined by means of multiplicative constant, spann also an algebra 
su{2). 

The spherical coordinates around the so axis are 

S\ = COS 02 COS 81, S2 = cos 6*2 sin0i, So = sin(#2)- 
We obtain a new pair of operators, that written in terms of the original variables (<j>i, <fo) are 

C^ o Zi ^ % = ± (cos (pi tan (fadfa — sm4>id^ 2 ) + 1/2) cosec<^>itan02 + (Z2+I/2) sin^icot^- 
These operators intertwin the Hamiltonians in the following way 

The 'global' operators (again 1/2 times the 'old' ones) close the third su(2). 

All these transformations {A ± , A, B ± ,B, C ± , C} (C = B - A) spann an algebra su(3). 
The s«(3)-Casimir operator is given by 

C = A + A- + B + B- + C + C~ + -A(A - 3/2) + -B(B - 3/2) + -C(C - 3/2). 

3 3 3 

We obtain u(3) by adding the central diagonal operator D := Iq — l\ — I2. 

The global operator convention can be adopted for the Hamiltonians H in the u(3)- 
hierarchy by defining its action on the eigenfunctions &(i 1: i2,i 3 ) of Hn lt i 2} i 3 -\ by 

Then, the Hamiltonian H can be expressed in terms of both Casimir operators, 

1 1 ^ 

H = AC--D 2 + —. (6.1) 

3 4 v ' 

Hence, the Hamiltonian can be written as a certain quadratic function of the operators 
{A^^B^jC^} generalizing the usual factorization for ID systems. 

In this way we have built an algebra u(3) of intertwining operators that, once fixed the 
initial Hamiltonian with parameter values (lo,h,h), generate a two-parameter Hamiltonian 
hierarchy 

{H(lo+m,h+m-n,l 2 +n)}} m,n G Z, 

where the points (Iq + m, l\ + m — n, I2 + n) lie on a certain plane D = d$. 



12 



One can prove that the eigenstates of this Hamiltonian hierachy are connected to the 
IUR's of u(3). Fundamental states $ annihilated by A" and C~ (simple roots of su(3)) 

Aj* e = Cj$ t = 

only exist when Zi = 0, 

$t(<l>i,<h) = N cos io+1 / 2 (</>i)sin 1 / 2 (0 1 ) C os io+1 (</> 2 )sin i2+1 / 2 ( ( /> 2 ), 
whit N a normalizing constant. The diagonal operators act on them as 

A$ t = -lo/2$ e , l = m, h = 0, m = 0,1,2,... 

(6.2) 

C$ t = -l2/2$ t , l 2 = n, n = 0,1,2,... 

This shows that is the lowest state of the IUR j± = m/2 of the subalgebra su(2) generated 
by {A^, A}, and of the IUR j2 = n/2 of the subalgebra su(2) closed by {C , C}. Such a 
su(3)-representation will be denoted (m,n), m,n G Z-°. The points (labelling the states) 
of this representation obtained from 3>£ lie on the plane D = m — n inside the ^-parameter 
space. 

The energy for the states of the IUR determined by the lowest state 1)6.2)1 with parameters 
(Z ,0,Z 2 ), is given <JSHJ) by 

E = (Z + h + 3/2) (Zq + l 2 + 5/2) = (m + n + 3/2)(m + n + 5/2). (6.3) 

Note that the IUR's labelled by (m, n) with the same value m + n are associated to states 
with the same energy. We call such IUR's a iso-energy series. This degeneration will be 
broken using the algebra so(6). 

6.2 The so(6)-hierarchy 

Making use of some relevant discrete symmetries, following the procedure of section [5~2*l the 
dynamical algebra u(3) can be enlarged to so(6). 

The Hamiltonian Hn 0l i lt i 2 ) 1)5.1)1 i s invariant under reflections in the parameter space 
(h, h,h) 

Jo : {lo,h,h) — > (—loih,h), h '■ (lo,hih) — * Go> — h,h), h '■ (lo,h,h) - > {loih,—fa) 

These symmetries, can be directly implemented in the eigenfunction space, giving by con- 
jugation another set of intertwining operators (jX = IiX Jj, i = 0, 1, 2) closing an isomorphic 
Lie algebra ju(3). They are (now labelled with a tilde) 

{A ± ,B ± , C±} C*}, 
For instance, the sets {^l 1 * 1 ,^!} and {Jr 1 ,^} close the two commuting su(2) of section l5~2l 
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The explicit expression for these new intertwining operators can be easily obtained in the 
same way as was done in ()5.15j) . They close a Lie algebra of rank 3, so(6). Instead of the 
six non-independent generators A, A, B, B, C, C it is enough to consider three independent 
diagonal operators Lo,Li,L,2 defined by Li^n at i u i 2 \ := k^(i Q i lt i 2 )' The Hamiltonian can be 
expressed in terms of the so(6)-Casimir operator by means of the 'symmetrization' of the 
it(3)-Hamiltonian ()6.1|) 

H = {A + ,A~}+{B + , B-}+{C + , C-}+{i+ A-}+{B + ,B-}+{C + ,C-}+L 2 +L 1 2 +L 2 2 +^. 

The intertwining generators of so(6) give rise to larger 3D Hamiltonian hierarchies 

i^('o+m+p,ii+m-n-p,« 2 +™)}' m,n,p £ Z, 

each one including a class of the previous ones coming from u(3). The eigenstates of these 
so(6)-hierarchies can be classified in terms of so(6) representations. 
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